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Extending the Utility of Partially Quenched QCD 
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It has been proposed that partially quenched chiral perturbation theory can be used together with partially 
quenched lattice QCD to determine the low-energy constants of real QCD. We point out a complication of this 
approach and show how it can be resolved. Partial quenching changes the chiral symmetry group and introduces 
a new independent four-derivative operator. We explore the effects of this new operator. We describe tree-level 
scattering processes to which it contributes, and show how they can be used, in principle, to deterimine its 
unknown coefficient. We also calculate its lowest-order (one-loop) contribution to charged meson masses and 
decay constants. In addition, we present the form of the analytic next-to-next-to-leading order meson mass and 
decay constant corrections, which are needed for lattice fits. 

Introduction and Overview 

Chiral perturbation theory (xPT) allows ana- 
lytic calculation of low-energy QCD processes in- 
volving the light pseudoscalar mesons. Results 
are given in terms of a number of undetermined 
coefficients, such as the Gasser-Leutwyler coeffi- 
cients and /,r, which we would like to determine 
from first principles. Lattice QCD provides a 
method for doing this as long as simulations are 
done at low enough quark masses that xPT (typ- 
ically at next-to-leading order) is a good approx- 
imation. In practice, however, lattice simulations 
are not done at the actual QCD values of the 
quarks masses because simulating light, dynam- 
ical quarks is computationally expensive. Thus 
there is a limited range of quark masses where 
both lattice simulations are feasible and xPT de- 
scribes the physics. 

This situation can be improved using the par- 
tially quenched (PQ) approximation, in which 
valence quarks (those which appear in external 
states) and sea quarks (those which appear in 
dynamical loops) are allowed to have different 
masses. Allowing valence and sea quarks to have 
different masses expands the parameter space 
available for determining the QCD low-energy 
constants, as shown in Fig. ^ 

It has been argued in Refs. (Il2) that although 
PQ QCD is unphysical, the properties of pseudo- 
Goldstone bosons (PGBs) can be expressed in 
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Figure 1. Partially quenched parameter space 

terms of the physical coefficients of the QCD 
chiral Lagrangian through next-to-leading order 
(NLO) in xPT. While the conclusions about NLO 
masses and decay constants are correct, the argu- 
ment is incomplete. We point out an overlooked 
subtlety in using PQ simulations to extract physi- 
cal QCD constants, and discuss the consequences. 

New Four-Derivative Operator in PQxPT 

Partial quenching introduces unphysical 
bosonic ghost quarks and therefore changes the 
low-energy symmetry structure The chiral 
symmetry group becomes a graded group con- 
taining both commutation and anti-commutation 
relations: 
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SU {Nvalence + N Sea\Nv alence) (1) 
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We observe that, because of the graded group 
structure, the four-derivative operator: 



Str(9^f79^[/t9^[/9''C/t) 



(2) 



is hnearly- independent, while in SU{N), N < 3, 
it is not. For comparison with unquenched QCD, 
it is convenient to use a hnear combination of 
operators that vanish in the unquenched SU (3) 
sector of the PQ theory, with an undetermined 
coefficient, Lpq: 

OPQ = LpQ{Str{d^Ud,U^d''Ud''U^) 
-istr(9^C/a^C/t)2 

+2Str(a^C/9^[/^9^C/9'^C/t)} (3) 
In summary, we find that this new operator: 

• contributes at NLO to "unphysical" scat- 
tering processes that cannot occur in un- 
quenched SU{3). 

• contributes at next-to-next-to leading order 
(NNLO) to PGB masses and decay con- 
stants. 

Because O-pQ vanishes in the unquenched 
SU{3) limit, it cannot contribute at tree-level to 
"physical" scattering processes that can occur in 
the unquenched SU{3) sector.^ For example, the 
process shown in Fig. [5] is an SU{3) process be- 
cause it involves only two quark flavors. The vari- 
ous quark line diagrams cancel each other out be- 
cause of relative minus signs and numbers of Wick 
contractions. However, O-pQ does contribute to 
the scattering process in Fig. ^ because the pro- 
cess involves four quark flavors. In the case of 
PQ QCD with Nsea = 3, these could be four 
valence quarks, three sea quarks and a valence 
quark, or other possibilities. Figure 0] also re- 
ceives contributions from OpQ because this pro- 
cess is unique to the PQ theory, involving both 
valence and ghost quarks. These examples show 
that, physically, the new particle content of the 

^Here we do not consider processes involving neutral 
mesons because double poles in the meson propagators 
make it unclear how to amputate diagrams and thus de- 
fine scattering amplitudes. 



PQ theory allows separation of new combinations 
of quark contractions, necessitating an additional 
operator in the Lagrangian. 
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Figure 2. "Physical" scattering process because 
it involves only two quarks 
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Figure 3. "Unphysical" scattering process be- 
cause it involves more than three quarks 





Figure 4. "Unphysical" scattering process be- 
cause it involves both valence and ghost quarks 

Like /tt and the Gasser-Leutwyler coefficients, 
LpQ can be obtained by fitting to lattice ex- 
periments. To illustrate this we calculate me- 
son scattering processes that only receive con- 
tributions from double supertrace operators in 
PQxPT, such as that shown in Figure 0] The 
0{p^) part of the scattering amplitude is: 

fM^32Li{pi-p3)ip2-p4) 

+ I6L2 (pi ■P2){P3 - Pi) + {Pl ■P4){P2 -Ps) 

+16LpQ {pi ■ P2){P3 ■ Pa) + {pi ■ P3){P2 ■ Pa) 
+ {Pi ■ Pa){p2 ■ P3) (4) 
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Li, L2, and LpQ can, in principle, be separately 
determined by combining this with other ampli- 
tudes and fitting to lattice data^. 

Because O-pg is a four-derivative operator, it 
does not contribute to masses and decay con- 
stants until NNLO order in PQxPT, and there- 
fore does not affect the conclusions of Refs. |ll2j 
As a first step in extending the program of 
Refs. |1I2| to NNLO, we calculate the lowest-order 
contribution of O-pQ to charged meson masses 
and decay constants 5'. To simplify calculations 
we consider only two valence quarks, A and B, 
and N sea quarks of equal mass. The charged 
meson, ttab, receives NNLO corrections from the 
new operator to its mass from Figure [S] and to its 
decay constant from Figure All of the meson 
diagrams receive contributions from the classes 
of quark line diagrams shown in Figure which 
illustrate the flow of quarks within the mesons. 
Quark line diagrams show the effect of partial 
quenching — valence quarks appear only as ex- 
ternal states because ghost quarks cancel them 
in loops. Quark "hairpins" remove the flavor sin- 
glet propagator to enforce the fact that SU {N\M) 
generators are traceless. Both the mass and decay 
constant corrections from the new operator van- 
ish in the unquenched SU{3) limit (mA = mB — 
nT-Sea, Nsea = 3) as they should. The expressions 
will appear in Ref. 0. 

Analytic NNLO Mass and Decay Constant 
Corrections 

Fits of present PQ lattice data require NNLO 
terms |6I7| . Full non-analytic NNLO calculations 
in PQxPT are not available; however, as an in- 
termediate step we have determined the form of 
the analytic NNLO terms. Although over a dozen 
operators in Cq contribute, only four linear com- 
binations of quark masses appear in the tree-level 
correction to the mass of ttab'- 
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+ asixA + XsJ +a'i(xA-XBj (5) 



Here are linear combinations of £g coefficients. 
Corrections to Jab have the same form with dif- 
ferent coefficients. Fits to PQ lattice data with 
^ms < msea < 1*^5 ^re consistent with this 
NNLO formula (7j. 
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Figure 5. PGB mass renormalization from O-pQ 
vv, vs, ss vv, vs, ss 
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Figure 6. PGB decay constant renormalization 
from O-pQ 





^Problems due to the lack of unitarity in the PQ theory 
HI do not appear until one-loop level, i.e. NNLO in xPT. 
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Figure 7. Quark line contributions to mass and 
decay constant renormalization 
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